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Abstract 

We consider a network of telescopes capable of scanning all the observable 
r-| ! sky each night and targeting Near-Earth objects (NEOs) in the size range 

fp-f of the Tunguska-like asteroids, from 160 m down to 10 m. We measure the 

performance of this telescope network in terms of the time needed to discover 
at least 50% of the impactors in the considered population with a warning 
time large enough to undertake proper mitigation actions. The warning times 
are described by a trimodal distribution and the telescope network has a 50% 
probability of discovering an impactor of the Tunguska class with at least one 
week of advance already in the first 10 yr of operations of the survey. These 
results suggest that the studied survey would be a significant addition to the 
current NEO discovery efforts. 
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1. Introduction 

In the design of a NEO survey there is a possible trade-off between cover- 
in g less sky to a deeper magnitude or more sky to a fainter one, as described 



by iTonryl (l201ll ) . In the present paper we denote the first strategy as "Deep 



Survey", and the secon d one "Wide Sur vey" . In the literature, the basic idea 



of a Deep S urvey is in iMorrisonI (Il992l) . while the idea of a Wide Survey is 



described in iHills and LeonardI ( 119951 ) . The choice between the two observing 
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strate gies is driven by the goals of the survey (IStokes et al.l . 120021 ) . Accord- 
ing to iMorrisoru ( 119921 ) , Deep Surveys such as the present American ones are 
more effective in reaching the completeness of the NEO population as they 
scan larger volumes of th e Near Earth space fo r a fixed absolute magnitude. 
As a matter of fact, iMainzer et al.l (120111 ) claim that more than 90% 



of objects larger than 1 km (first Spaceguard goal) have been discovered 
so far, predominantly by US surveys. Nevertheless, these surveys are not 
optimized for detecting imminent, relatively small impactors from 10 to 160 
m diameter, which may s till cause very impo rtant damages and losses on 
the ground. As shown by iBrown et al.l (120021 . Fig. 4), th e energy r e leased 
by such impactors ranges from 50 to 10^ kT and already iMorrisonI (119921 ) 
discusses the substantial local damage that a Tunguska-sized impactor can 
inflict to a populated area. 

The reason why deep surveys are not suited for imminent impactors is 
that their observing strategy is to cover the same area in the sky after a 
few days, and to take only a minimum number of images. This impairs 
the successful identi fication of o b jects that are going to impact within a few 
days. For instance, IVeres et al.l (120091 ) prove that Pan-STARRS would not 
have been able to collect enough detections to compute an orbit for 2OO8TC3, 
a ~ 5 m asteroid which impacted on Earth on October 2008. 

To deal with imminent impacts, a more effecti ve st rategy is the Wide 
Survey, as demonstrated by lHills and Leonard! ( 119951 ) and lTonryl (120 111 ). This 
kind of survey provides a more responsive NEO impact warning system and 
thus could nicely complement the current NEO discovery and cataloging 
strategy of the US programs. 

In the present paper we measure the performance of an assumed Wide 
Survey design through a simulation of a 100 yr time span of operations. We 
deal with small impactors, i.e., those with absolute magnitude between 22 
and 28, and measure the time it takes to reach a 50% threshold for the frac- 
tion of objects discovered with warning time sufficient to undertake proper 
mitigation actions. 



2. Blind time 

An impactor can arrive from almost anywhere in the sky. In particular, 
if the impactor comes from the direction of the Sun, it will be most likely 
not detectable in the last days before its fall. This implies that such an 
object should be discovered at a previous apparition, if at all possible. Thus, 
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we are led to the possibility that, after the beginning of the operations of a 
survey, there is no chance for the potential impactor to be discovered before 
its fall; i.e., the survey is "blind" for this specific impactor. Such a situation 
depends on several factors, including the orbit and absolute magnitude of the 
impactor, and the parameters characterizing the survey (limiting magnitude, 
sky coverage, cadence, etc.). 

Given an impactor, we define the "lead time" as the interval of time 
between the first orbit determination and the time of impact. According to 
the size of the impactor, the lead time should be large enough to undertake 
the required mitigation actions, i.e., the larger impactor the more time is 
necessary for the mitigation. In absence of specific information about the 
albedo and the shape of an imminent impactor, the size can be inferred from 
its absolute magnitude H. Thus, we define the minimum required lead time 
as a function of if, using the following constraints: 

• a minimum lead time of 30 days for objects of H = 22; 

• a minimum lead time of about one week for Tunguska-sized impactors 
{H = 24.5). 

We adopt the following function 

t(d) = cie-"^(^-22) , ci = 30 d , C2 = 0.5 . (1) 

that fulfills the above constraints. Figure [1] shows t as a function of H. 

It is important to point out that this simple law is tailored to the pop- 
ulation used in our simulations, with H ranging between 22 and 28 (see 
Sec. 13. 2p . For these objects, which are the target of the Wide Survey de- 
scribed in this paper, the mitigation actions to be undertaken are essentially 
orbit improvement and evacuation. Dealing with bigger objects requires a 
different approach and mitigation strategy, thus Eq. ([T]) should be replaced 
with a different model, possibly involving a larger number of parameters. 

Given an NEO population, we define as "blind time" for a given survey 
and a given absolute magnitude H, the time between the start of operation of 
the survey and the moment at which 50% of the impactors with magnitude H 
have a lead time larger than the minimum threshold defined by Eq. ([T]). The 
blind time can be used as an indicator of the performances of a given survey. 
As a metric, the blind time is a variant of the time required for a survey 
to discover 90% of a defined population. The usefulness of this definition is 
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Figure 1: Minimum lead time as a function of the absolute magnitude H . 

that, when deahng with small but numerous NEOs, the time scale for a 90% 
completion is very long and uncertain, due to the poor modeling of the small 
object population. 

3. The simulation 

Hereafter we describe our assumptions on the optical network, on the 
impactor population, and the orbit determination process used in the simu- 
lation. 

3.1. Optical network 

For the optical sensors we assume the use of the innovative fly-eye tele- 
scope design, having the following main characteristics: 

• an equivalent aperture of 1 m; 
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• a FoV of 45 deg^ (6.66° x 6.66°); 

• high efficiency CCDs (80-90%) with very fast read-out times (~ 2s) 
and very good cosmetics (~ 99%); 

• a fill-factor ~ 1, that is the ratio of the effectively detected area of the 
FoV and the FoV; 

• a minimum elevation of 15° above the horizon. 

The above assumptions on the sensor hardware require a significant effort 
in both technological developme nt and resources. T he concept design of the 



assumed telescope is described in lCibin et al.l (120111 ). However, we are aware 



that these assumptions need to be validated by further studies when the 
telescope is actually available, at least in a prototype phase. A discussion on 
this is beyond the scope of this paper. 
For the telescope network we assume: 

• One equivalent dedicated survey telescope in the northern and one in 
the southern hemisphere. 

• The northern telescope covers the northern hemisphere of the celestial 
sphere, while the southern telescope covers the southern hemisphere. 

• One dedicated follow-up telescope in the northern and one in the south- 
ern hemisphere, typically 30° West of the survey telescopes. 

• The images are processed locally in real time, included the astrometric 
reduction, and the data are made available to the scientific community 
in less than two hours. Therefore, the dedicated follow-up telescopes 
can be triggered to follow the newly discovered objects. 

With these assumptions each equivalent telescope can take about 766 images 
for an average 10 hour night. This corresponds to a total of about 34450 
deg^ which is equivalent to 17225 deg^ of the celestial sphere taken twice per 
night. 

For the observing strategy we assume: 

• Observations that cover 36400 deg^ (^ 88% of the celestial sphere), 
corresponding to all the visible sky except the regions with solar elon- 
gation less than 40°. 
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• The regions of the sky within 30° of the Moon or within 15° of the 
galactic plane are not covered by the telescopes due to the increase of 
the sky background. Therefore, the effective visible sky ranges between 
22987 deg^ ~ 56% of the celestial sphere (when the forbidden regions 
around the Sun, the Moon and the galactic plane do not overlap) and 
34348 deg^ ~ 83% (when the intersection between the forbidden regions 
is maximized). On average, each telescope covers between 11500 and 
17200 deg\ 

• A limiting magnitude Vum = 21.5, corresponding to ^ 45 s of exposure 
time, for the survey mode, and Vum = 23 for the follow-up mode. 

• Coverage of the visible sky at least two times per night. 

In a real system the number of telescopes has to be increased and could be 
between 5 and 6. Indeed, to deal with the cloud coverage and meteorological 
correlations we need a minimum of two survey telescopes per hemisphere 
widely spread in longitude. Furthermore, to increase the detection efficiency, 
a higher number of detections may be necessary and also this can be achieved 
with a higher number of survey telescopes. 

3.2. Impactor population 

In our simulation we use t h e pop ulation of 4950 synthetic impactors de- 



scribed by IChesley and Spahrl (120041 ) which have impacts in a time frame of 



100 yr starting from July 2009. This impacto r population is selected within 



the population model by iBottke et al.l (120021 ). Figure [2] shows the distribu- 



tion of semimajor axis, eccentricity and inclination. The majority (68%) of 
the objects has a perihelion between 0.8 and 1 AU. 

We assign a fixed value of the absolute magnitude H to all the asteroids 
and we repeat the simulation for integer values of H ranging between 22 and 
28, roughly corresponding to diameters between 160 and 10 m. We choose 
this simulation strategy to measure the performance of the proposed network 
as function of the size range of the asteroids. 

To obtain the sky distribution of the impactors shortly before the event. 



it is useful to plot their radiants (see Appendix A ). Figure |3] shows the dis- 
tribution of the radiants for 4 465 Earth impactors in such a represent at ior0. 



The analytical procedure to compute the radiants assumes a circular orbit for the 
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The sky distribution of impactor radiants is far from uniform and is a con- 
sequence of the a-e-i distribution of the impactor population. The fraction 
of radiants with a solar elongation larger than 40° — the minimum elonga- 
tion from the Sun at which the assumed survey can observe — is 80.1% of the 
whole sample. It is worth noticing that, for the impactors with radiant within 
40° of the Sun, a detection at an apparition before the one corresponding to 
the impact is the only chance to have a long lead time. 

3.3. Methodology 

We split the impactor population in 10 bins according to the impact epoch 
with respect to the beginning of the simulation. For example, the first bin 
contains objects impacting within 10 yr, the second one objects impacting 
between 10 and 20 yr, and so on. Such a binning allows us to measure the 
performance as a function of the time from the start of the survey. 

For each object we generate a list of observations according to the assumed 
configuration, the performance of the optical network, and the visibility con- 
straints. The simulation provides one tracklet (see later) per night for survey 
telescopes and up to two tracklets per night for follow-up telescopes. The 
follow-up observations are triggered once the object has been detected by 



Earth. 485 among the objects of Cheslev and Spahil ( 2004 ) have either a(l — e) > 1 or 



a(l + e) > 1, so that they are excluded from the analytical computation. 
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Figure 3: The radiant distribution of the simulated impactors o f ICheslev and Soahil ()2004l) . 
The radiants are shown in an equal area projection of the sky centered on the opposition; 
the angular coordinates are ecliptic longitude minus the longitude of the Sun, and ecliptic 
latitude. The bold lines refer to 40° of solar elongation. 



the survey telescope(s), with a minimum delay time of two hours. For each 
simulated observation Gaussian noise with a standard deviation of 0.3 arcsec 
is added to the astrometric position. Similarly a Gaussian noise is added to 
the magnitude estimate for the detection, but in this case the noise is split 
in a correlated component (for the light curve effect) of 0.2 magnitudes and 
a random component, again of 0.2 magnitudes. 

The tracklet is the atomic information for a moving asteroid, consisting 
of a small number (2-5) of detections in different images of the same field, 
taken at moderately short intervals of time (15 min to 2 hours). A tracklet 
normally provides an amount of information which can be described by 4 
scalar quantities (two ang les and two angular rates) , therefore such detections 



do not imply discovery ( Milani et al. . 2007f ). We consider as discovered a 



moving object belonging to the solar system only when enough information 
has been accumulated to establish its dynamical properties, that is by means 
of a heliocentric orbit, for which at least 6 scalar quantities are required. 

The orbit determination process starts by selecting n-tuples of tracklets 
which could belong to the same object. Then, for each of them a prelimi- 
nary orbit c ompatible wit h all t he tr acklets is computed, using t he methods 
described in iMilani et al.l (12004[ 1 and iMilani and Gronchil ( 120101 . Chapters 7 
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and 8). Thereafter, the prehminary orbit is used as first guess in a differ- 
ential correction procedure, which usually converges to a least squares fit 
orbit. If the orbital fit satisfies suitable quality control conditions, this can 
be considered a real object. 

To perform the simulation we set up a data center architecture, by in- 
gesting observational data day by day. Each time new observations become 
available we update the previously known orbits and compute the new ones, 
corresponding to newly discovered objects. 

It is important to point out that a more realistic simulation should take 
into account the presence of Main Belt background asteroids, which increases 
the computational load and the rate of occurrence of false identifications. 
However, observations of known Main B elt asteroids should be filtered before 
looking for new objects, as discussed in iMilani et al.l (120 121 ). 

For completeness, our procedure could include the risk assessment for 
each simulated impactor including the explicit computation of an impact 
probability. Such a procedure would requir e a very large amount of CPU 
time and was performed in a previous study (IFarnocchia et al.l . l201l[ ) using a 
subsample of the present population over the first 20 yr of survey operations. 
That paper shows that, in more than 99.5% of the cases, the availability of an 
orbit involving at least 4 tracklets for a newly discovered potential impactor 
is accompanied by the successful cor nputation of an ira pact probability by 
using the CL0M0N2 software robot (IMilani et al.l . 120051 ). Such a percentage 
is so high as to make not cost effective the computation of the entire impact 
monitoring chain in the present case. Thus, we stipulate that an impactor is 
considered as discovered when an orbit from at least 4 tracklets is computed. 



4. Results and Discussion 

The main outcome of the simulation is shown in Fig. HI The left panel 
shows the differential discovery completeness as a function of time from the 
survey beginning and absolute magnitude. For each 10 yr bin, the differ- 
ential discovery completeness is defined as the ratio between the number of 
impactors discovered and those impacting in that time frame. For example, 
the differential discovery completeness for H = 24 and the bin from 20 to 
30 yr is given by the fraction of objects discovered before the impact among 
those impacting between 20 and 30 yr with respect to the beginning of the 
survey operations. 
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Figure 4: Left: differential discovery completeness as function of the impact date and 
absolute magnitude. Right: same as left panel, considering successful a discovery with 
lead time greater than the minimum threshold defined by Eq. ([T]). 

For impactors with H = 23 a 70% discovery completeness is achieved in 
the first decade, while the 90% threshold is reached after two more decades. 
For impactors with H = 25 a 60% discovery completeness is achieved after 
two decades, while the 90% threshold is almost reached at the end of the 
simulation. For H = 28 the discovery completeness starts slightly below 
40%, and increases slowly during the next decades as expected, due to the 
small size of the objects. Notice that the completeness is greater than 50% 
already from the start for H = 26. 

The right panel shows the differential completeness for a lead time given 
by Eq. ([1]). This means that a discovery is considered successful only if it 
takes place sufficiently early, allowing the mitigation actions appropriate for 
the size of the impactor. For impactors with if = 23 a ~ 60% completeness 
is achieved in the first decade, while the 90% threshold is reached after 30 yr. 
For if = 25 a ^ 50% discovery completeness is achieved in the second decade, 
while the 90% threshold is not yet reached at the end of the simulation. For 
smaller asteroids the discovery completeness is larger than 25% from the 
beginning of the operations and overcomes the 50% threshold after more 
than 70 yr. 

We can sum up the results in terms of blind time, that is the intersection 
between the 50% completeness and the H = const curves in Fig. |H The 
survey simulated in this paper would have a blind time of about 20 yr for 
imminent impactors with if = 25. Note that a Tunguska-sized (if ~ 24.5) 
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object impacting within 10 yr from the start of the survey would have a 
> 60% probabihty of being discovered and would have a lead time larger 
than 1 week with a probability ~ 45%. For smaller impactors the blind time 
increases up to ~60 yr for H = 28. 




Figure 5: Histogram of the lead times for different values of H =22 (top left), 24 (top 
right), 26 (bottom left), 28 (bottom right). The vertical lines denote, from left to right, 1 
day, 1 week, 1 month, and 1 yr. Note that the scales are different. 

Figure |5] shows the distribution of the lead time for different values H. 
As expected, the lead time strongly depends on the value of the absolute 
magnitude. A clear trimodality is visible in all the panels: either the object 
impacts without being discovered (left bar), or is discovered during its last 
apparition (central peak), or at a previous apparition (right peak). Tabled] 
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details the fractions of objects in each peak for a fixed value of H. Most 
of the impactors with H = 22 are discovered during a previous apparition 
with respect to the impact. As H increases there are more and more cases of 
objects either discovered during the last apparition or not discovered at all. 




Years before impact 



Years before impact 



Figure 6: Left: integral discovery corapleteness as function of the impact date and absolute 
magnitude. Right: same as left panel, considering successful a discovery with lead time 
greater than the minimum threshold defined by Eq. p]). 

To conclude this discussion we report the integral completeness achieved 
by our simulated survey. The integral completeness is computed by the 
weighted sum 

Comp(i7 <H)= WH,Comp{Hi) / wh^ 

22<Hi<H 22<H,<H 

where Comp(if) is the completeness for a fixed absolute magnitude H. To 
get a more reahstic result we keep into account a power law distribution for 
the number of aste roids at a given absolu t e mag nitude. Consistently with 
Bottke et al.l (120021 ) and IStuart and Binzell (120041 ) , the computed efficiencies 
for a fixed absolute magnitude value H are given the weights: 

The results are summarized in Fig. |6l 



5. Conclusions 

We simulated the operations, over a time span of 100 yr, of a Wide Survey 
capable of covering all the sky at solar elongation larger than 40°, down to 
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H 


undiscovered 


last apparition 


previous apparition 


22 


4.2% 


8.9% 


86.9% 


23 


8.1% 


13.6% 


78.4% 


24 


15.0% 


18.5% 


66.5% 


25 


22.3% 


25.8% 


51.9% 


26 


31.4% 


31.0% 


37.6% 


27 


41.3% 


33.0% 


25.7% 


28 


52.1% 


31.1% 


16.8% 



Table 1: Percentages of impactors not discovered (2nd column), discovered during the last 
apparition (3rd column) or discovered at a previous apparition (4th column), as a function 
of (1st column). 

apparent magnitude 21.5, with a nightly cadence. The survey includes the 
operation of follow-up with a limiting magnitude 23.0. 

The goal of the simulation was to compute the "blind time" of the sur- 
vey, i.e., the time between the start of survey operations and the moment at 
which 50% of the impactors at a given magnitude are discovered and their 
orbits determined, with an advance large enough to allow undertaking the 
appropriate mitigation actions. In fact, our modeling allowed us to compute 
a realistic distribution of the lead time, i.e., the interval of time between the 
first orbit determination and the time of impact. The distribution shows 
a trimodality corresponding to 1) undiscovered objects, 2) objects discov- 
ered at the impact apparition, and 3) objects discovered during a previous 
apparition. 

The survey discussed in this paper can efficiently deal with Tunguska- 
sized impactors for which the blind time is about 10 yr. This means that, 
already in the first ten years of survey operations, there would be 50% prob- 
ability of discovering such an impactor at least one week before impact. The 
pre-impact discovery of Tunguska- sized impactors was not an original goal 
of past NEO surveys but the case of 2008TC3 and the present paper show 
that it is becoming within reach. 

In a future paper, we plan to compare the performances of the discussed 
Wide Survey and a state of the art Deep Survey thus allowing us to quan- 
titatively measure the contribution of the Wide Survey to the current NEO 
search programs. 
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Appendix A. Impactor radiants 



Following IValsecchi et al.l (119991 ) it is possibl e to e s tablis h that in the 
framework of Opik's theory of close encounters (lOpikl . Il976l ) the angles 9 
and (f) define the direction opposite to that from which an Earth impactor 
seems to arrive (the so-called radiants in meteor astronomy). These angles 
can be computed from the orbital elements a, e, i of the impactor (a must be 
given in units of the orbital semimajor axis of the planet, for Earth impactors 
in AU) as follows: 



cos 9 



sin 9 



smi 



COS( 



COS? 



1 



3 — 1/a — 2^^0(1 — e^) cosi 



^2- l/a-a(l -e2) 



cos^ i 



3 — 1/a — 2^J a{l — e^) cosi 

v/2- l/a-a(l -e2) 
a/2 — 1/a — a(l — e^) cos^i 

2 — 1/a — a(l — e^) cos^ i ' 



where in the expression for sin0 the upper sign applies to collisions in the 
post-perihelion branch of the orbit, and in that for cos0 to collisions at the 
ascending node. 

Thus, to each triple a, e, i correspond four encounter geometries, all 
characterized by the same value of 6', that differ for the quadrant of 0; this, 
in turn, can be computed from the orbital elements of the impactor, as shown 
in Table \K?2\ 

Meteor radiants are often plotted using an equal area projection centered 
on the direction of the Earth motion, with the ecliptic as reference plane; 
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sin q 


^ > 


cos</ 


) > 


u 


+ f = 





< / < TT 


sin q 


;> < 


COS 


i > 


UJ 


+ / = 





71 < f <2n 


sin q 


i > 


COS(/ 


) < 


u 


+ / = 


TT 


< / < TT 


sin q 


!) < 


COS q 


i < 


UJ 


+ / = 


71 


TT < / < 27r 



Table A. 2: The relationship between the quadrant of (p and the orbital elements ui and / 
of the impactor at collision. 

in this case, starting from 6 and computed from the meteoroid orbits, the 
radiant coordinates are simply given by A + vr and —(3 {/3 = defines the 
echptic plane, and A = /3 = is the direction of the Earth motion), with A 
and /3 given by 

sin ^ sin cos6' 

smp = smt^cos0 , smA = — , cosA = -. 

cos p cos p 

For an asteroidal Earth impactor, it is useful to use the same representation, 
but in this case centered on the opposition point, something that is obtained 
by rotating A by 7r/2 in the appropriate direction. 

Appendix A.l. An analytical expression 

The non-uniform sky distribution shown in Fig. [3] can be exploited in 
prioritizing the sky coverage for a survey aimed at detecting very close Earth 
approachers and impactors. 

In this respect, an important quantity to take into account is the angular 
distance of the radiant from the Sun, that we will denote with a, since with 
ground based optical telescopes it is practically impossible to observe at 
values of this quantity smaller than some practical limit. Thus, in order to 
evaluate the efficiency of an impactor-aimed sky survey, it is important to 
establish what fraction of impactor radiants lies below this angular distance 
from the Sun; to this purpose, it is nec essary to discuss the geometric setup 



of Opik's theory IValsecchi et al.l ( 119991 ). 



We use a reference frame centered on the Earth, with the 2;-axis perpen- 
dicular to the plane of the ecliptic, the y-axis in the direction of the Earth 
velocity and the x-axis pointing away from the Sun, which is located at 
X = —1, y = 0, z = 0; in this frame, the unperturbed geocentric encounter 
velocity U of the NEO has components 

{Ux, Uy, Uz) = {U sin 6 sin 0, U cos 6, U sin 6 cos (p) 
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where U is the magnitude of the velocity vector. With these definitions, the 
cosine of the angle between U and the x-axis is simply given by sin 6 sin 0; 
this angle, in turn, is equal to 7r — a since, as noted before, 9 and 4> define 
the direction opposite to the radiant. 

Thus, if amin is the minimum angular distance from the Sun that the 
survey can reach, the radiants that are not observable by it are characterized 

by 

cos amin < sin ^ sin 0. 

Since for each triple a, e, i we have four possible associated radiants, two of 
which with sin0 > and the other two with sin0 < 0, the consequence of 
the above inequality is that for o"mm < 7r/2 we have the following two cases: 

• for I cos amin I > I sin 6' sin 0| , all four radiants associated to a given triple 
a, e, i are observable; 

• for I cos CTmin I < | sin^sin0|, the two radiants for which sin0 < are 
observable, while the other two are not. 

Given a population of impactors, it is then possible to compute the fraction 
F of radiants that have a > amin] if h is the fraction of the population 
characterized by | cos a min I ^ I sin sin (j) I , then 

„ , l-h l+h 
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